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ABSTRACT 

We report the first clear evidence for the simultaneous presence of a low frequency break and a 
QPO in the fluctuation power spectrum of a well known ultraluminous X-ray source (ULX) in M82 
using long XMM-Newton observations. The break occurs at a frequency of 34.2tg mHz. The QPO 
has a centroid at vqpo — 114.3 ± 1.5 mHz, a coherence Q = vqpo / ^ v fwhm — 3.5 and an amplitude 
(rms) of 19% in the 2 — 10 keV band. The power spectrum is approximately flat below the break 
frequency and then falls off above the break frequency as a power law with the QPO superimposed. 
This form of the power spectrum is characteristic of the Galactic X-ray binaries (XRBs) in their high 
or intermediate states. M82 X-l was likely in an intermediate state during the observation. The EPIC 
PN spectrum is well described by a model comprising an absorbed power-law (L ~ 2) and an iron 
line at ~ 6.6 keV with a width a ~ 0.2 keV and an equivalent width of ~ 180 eV. Using the well 
established correlations between the power and energy spectral parameters for XRBs, we estimate a 
black hole mass for M 82 X-l in the range of ~ 25 — 520MQ including systematic errors that arise 
due to the uncertainty in the calibration of the photon spectral index versus QPO frequency relation. 
Subject headings: accretion, accretion disks — stars: individual (M82 X-l) — X-rays: stars 



1. INTRODUCTION 

Ultra-luminous X-ray sources (ULXs) are extra- 
nuclear point X-ray sources with luminosities exceeding 
~ 10 39 erg s _1 . These objects were first discovered with 
Einstein observations (reviewed by Fabbiano 1989). The 
nature of ULXs continues to be an enigma, since their 
adopted isotropic high energy output surpasses the Ed- 
dington limit of even the most massive stellar mass black 
holes (BHs), sometimes by large factors. Several models 
have been proposed to explain the high luminosities of 
ULXs. The most popular is the "intermediate mass black 
hole (IMBH)" with mass M BH ^ 10 2 - 10 4 M Q (e.g., 
Colbert & Mushotzky 1999, hereafter CM99) bridging 
the gap between stellar mass BHs in XRBs and super- 
massive BHs in active galactic nuclei 3 . Other popular 
models include XRBs with anisotropic emission (King 
et al. 2001), beamed XRBs with relativistic jets di- 
rectly pointing towards us i. e., scaled down versions 
of blazars (Mirabel & Rodriguez 1999), and XRBs with 
super-Eddington accretion rates (Begelman 2002). Sev- 
eral observations suggest that ULXs may be similar or 
scaled up versions of XRBs. Discovery of orbital mod- 
ulations from several ULXs (Bauer et al. 2001; Sugiho 
et al. 2001) implies their binary nature. Recent Chan- 
dra and XMM-Newton observations of ULXs show soft 
X-ray excess emission which has been interpreted as the 
emission from accretion disks with temperatures in the 
range - 100-500 eV (see Shrader & Titarchuk 2003 and 
review by Miller & Colbert 2003 and references therein) . 
Observations of spectral transitions between low/hard 
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and high/soft states in two ULXs in IC 342 (Kubota et 
al. 2001) further demonstrate their similarity with the 
XRBs. The recent observation of a break at a frequency 
of ~ 2.5 mHz in the power density spectrum (PDS) of 
the ULX in NGC 5408 suggests a BH mass of - 100M© 
(Soria et al. 2004). 

The enigmatic nature of ULXs can be understood by 
determining their BH masses. The most reliable method 
is to measure the mass function through the secondary 
mass and orbital parameters such as velocities, period, 
orbit size etc., which can be measured only if the sec- 
ondary is optically identified. There are two cases of op- 
tical identification (Kuntz et al. 2005; Liu et al. 2004); 
however, the orbital parameters of any ULX binary sys- 
tem are yet to be measured. Good progress can still 
be made by establishing a direct physical connection be- 
tween ULXs and XRBs. This can be done by compar- 
ing the characteristic time scales such as those associ- 
ated with the low frequency quasi-periodic oscillations 
(QPOs). However, this comparison is ambiguous with- 
out the information about a low frequency break, the 
shape of the PDS and energy spectrum. 

M82 X-l is one of the brightest ULX (CXO M82 
J095550. 2+694047, source 7 in Matsumoto et al. 2001) 
that is well suited for a determination of the shape of 
the PDS. Chandra High Resolution Camera observations 
showed that the source is unresolved and off-nuclear 
(Matsumoto et al. 2001). The ULX has no optically 
bright counterpart (Kaaret et al. 2001): there is a likely 
IR counterpart which is a super star cluster (Kaaret et 
al. 2004). Recently, Strohmayer & Mushotzky (2003), 
hereafter SM03, discovered a QPO in M82 X-l at a fre- 
quency of 54 mHz based on a 27 ks XMM-Newton ob- 
servation. They also found QPOs in the 50 — 100 mHz 
frequency range in the RXTE data. Fiorito & Titarchuk 
(2004), hereafter FT04 apply a new method to determine 
the BH mass of M82 X-l. The method uses the spectral 
index-QPO low-frequency correlation that has been re- 
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cently established in BH XRBs GRS 1915+105, XTE 
J1550-564, 4U 1630-47, and others (Titarchuk & Fior- 
ito 2004, hereafter TF04) . Using scaling arguments and 
the correlation derived from Galactic BHs, they conclude 
that M82 X-l is an intermediate mass BH of order 1000 
Mq. We revisit the BH mass estimate in M82 X-l using 
new data on the power spectrum and energy spectra. 

2. OBSERVATION, DATA REDUCTION & ANALYSIS 

M82 was observed using XMM-Newton on 2001 May 
5 for 30 ks and again on 2004 April 21 for 103 ks. The 
first observation led to the discovery of a 54 mHz QPO 
from the ULX (SM03). Here we consider the EPIC data 
from the second observation only. The EPIC PN and 
MOS cameras were operated in the full frame mode us- 
ing the medium filter. We used the SAS version 6.1 and 
the most recent calibration data base to process and fil- 
ter the event data. Examination of the background rate 
above 10 keV showed that the observation is completely 
swamped by the particle background after an elapsed 
time of 71.5 ks and this latter period was therefore ex- 
cluded from the rest of the analysis. 

In the central regions of M82, Chandra images revealed 
diffuse emission and several point sources, M82 X-l be- 
ing the brightest among them (Matsumoto et al. 2001). 
The diffuse emission is evident in the EPIC images but 
the point sources are not resolved. Following SM03, we 
extracted PN and MOS events using a 18" circular re- 
gion around the bright point source whose position is 
consistent with the Chandra position of M82 X-l. Above 
2 keV, the diffuse emission contributes only < 10% to 
the point source inside the 18" circular region. At the 
faintest flux level of M82 X-l, nearby point sources con- 
tribute < 30% to the 0.5 — 10 keV X-ray emission in the 
18" region (Matsumoto et al. 2001; SM03). 

2.1. The power spectrum 

For temporal analysis, we combined the PN and MOS 
data and used the continuous exposure of 70.2 ks during 
which both the PN and MOS cameras operated simulta- 
neously. We calculated a power density spectrum (PDS) 
using the background corrected PN+MOS light curves 
sampled at 0.5 s. Figure ^ (left) shows the 2 — 10 keV 
PDS of M82 X-l rebinned by a factor of 1024 yielding 
a frequency resolution of 7.8 mHz. The PDS continuum 
is approximately flat at low frequencies below ~ 30 mHz 
and then falls off approximately following a power law 
up to ~ 200 mHz where the white noise arising from the 
poisson errors starts to dominate the PDS. There is a 
prominent QPO with its peak frequency near 114 mHz. 

We fitted the PDS with two models: (i) a broad 
Lorentzian for the continuum and a narrow Lorentzian 
for the QPO, and (ii) a broken power law (BPL) for 
the continuum and a Lorentzian (L) for the QPO. We 
also used a constant to account for the poisson noise. 
Both models provided statistically acceptable fits, giv- 
ing a minimum \ 2 °f 128.6 for 121 degrees of freedom 
(dof) and 127.9 for 120 dof for the double Lorentzian and 
BPL+L models, respectively. The errors on the best-fit 
PDS model parameters, quoted below, are at a la level. 
The double Lorentzian model resulted in a QPO centroid 
frequency vqpo = 114.6 ± 1.5 mHz, a width vfwhm = 
31.3±2 1 mHz, an amplitude A QPO = 0.037.4 ± 0.0024 
and the broad Lorentzian centroid vq = H.4^4'4 mHz, 



a width vfwhm — 60.2lg3 mHz and an amplitude 
A = 0.038 ± 0.0034. The BPL+L model resulted in 
a QPO centroid frequency vqpo = 114.3 ± 1.5 mHz, a 
width vfwhm = 32. 7~t 2 '_ 6 mHz, an amplitude Aqpo — 
0.038 ± 0.0024, and broken power-law indices of Ti — 
O.lllg'Jg below the break frequency v^ = 34.2^2'g mHz 
and T 2 = 2.32±gjf above the break frequency, and an 
amplitude Abpl = 0.36 ± 0.03 at a reference frequency 
of 10 mHz. The total integrated power (0.001 — 1 Hz) 
and the QPO power expressed as rms/mean are 23% 
and 19%, respectively. The use of a broken power law is 
an improvement (A\ 2 = —24.8 for two additional param- 
eters) over a simple power law at a statistical significance 
level of > 99.99% based on the maximum likelihood ra- 
tio test. Thus the break in the continuum of the power 
spectrum is real. The best-fit BPL+L model is shown in 
Fig-H^e/i) as a thick line. 

For comparison we have also calculated the PDS of an 
XRB XTE J1550-564 using the RXTE observation of 10 
September 1998. XTE J1550-564 shows a large range in 
its low frequency QPO centroid (vqpo ~ 0.08 — 18 Hz; 
Sobczak et al. 2000a). Our choice of the particular 
RXTE observation relies on the fact that XTE J1550- 
564 showed a power-law (PL) photon index of T ~ 2.0 
on 10 September 1999 (Sobczak et al. 2000b), which is 
very similar to the 3 — 10 keV photon index of M82 X-l 
(see below). Fig.\Hright) shows the PDS of XTE J1550- 
564 calculated from the PCA light curve sampled at 
0.125 s which was obtained from the RXTE public data 
archive. We used a model comprising a broken PL and 
two Lorentzians to fit the PDS. This model resulted 
in a statistically acceptable fit (minimum \ 2 — 94.2 
for 118 dof). The broken PL has best-fit parameters: 
Ti = 0.03 ± 0.006 below the break frequency Vf, — 
0.27 ± 0.007 Hz, T 2 = 0.89 ± 0.01 above the break and 
an amplitude Abpl — 0.026 ± 0.0005 at a reference fre- 
quency of 10 mHz. The two QPOs have centroid frequen- 
cies vqpo = 1.034 ± 0.004 and 2.055 ± 0.016 Hz, widths 
vfwhm = 0.125 ± 0.010 and 0.39 ± 0.05 Hz and ampli- 
tudes Aqpo = 0.026 ± 0.0005 and 0.0046 ± 0.0003 for 
the fundamental and first harmonic, respectively. The 
total integrated power (0.01 — 4 Hz) and the QPO (fun- 
damental) power expressed as rms/mean are 23% and 
16%, respectively. 

2.2. The energy spectrum 

We extracted PN and MOS spectra using a 18" cir- 
cular region centered at the position of M82 X-l. We 
also extracted PN and MOS background spectra using 
nearby circular regions free of sources. We created ap- 
propriate response files using the SAS tasks rmf gen and 
arfgen. The spectra were grouped to a minimum of 
20 counts per spectral channel and analyzed with XSPEC 
11.3. The errors on the best-fit spectral parameters are 
quoted at a 90% confidence level. An absorbed power 
law (PL) model fitted to the PN data showed strong soft 
excess emission below ~ 3 keV due to the presence of 
the diffuse soft X-ray emission in the source extraction 
region. We ignored the data below 3 keV and performed 
the spectral fitting of the PN and MOS data separately. 
A simple PL model poorly describes both the PN and 
MOS data. We show the ratio of the PN data and 
the best-fitting absorbed PL model in Figure Gla). A 
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prominent iron Ka line at ~ 6.6 keV is evident. Ad- 
dition of a Gaussian line (GL) to the PL model im- 
proved the fit significantly (A^ 2 = —71.7 for three addi- 
tional parameters) and resulted in a good fit (minimum 
X 2 = 898.9 for 944 dof. The Gaussian line has a centroid 
energy E Une = 6.61±° o8 keV > a width a = 278 -iil eV ; 
a line flux fu ne — 1.9^q'q x 10 -5 photons cm -2 s _1 and 
an equivalent width of ~ 180 eV. The best-fit pho- 
ton index is 2.00^q'qj and the the absorption column is 

N H = 2.5±£{j x 10 22 cm" 2 . The observed 2-10 keV 
flux is 9.6 x 10~ 12 erg cm~ 2 s^ 1 and the correspond- 
ing luminosity is 1.7 x 10 40 erg s _1 assuming a distance 
of 3.9 Mpc (Sakai & Madore 1999). In Figure EJd, we 
show the best-fit PL+GL model to the EPIC PN data. 
The PL+GL model fitted to the MOS data resulted in 
slightly flatter PL {T x = 1.80lg;J£), an absorption col- 
umn Ns = 1.8^0 3 x 10 22 cm~ 2 and a weaker iron line, 
EW ~ 94 eV. 

3. DISCUSSION 

Based on long XMM-Newton observations, we show 
clear evidence for the simultaneous presence of a low fre- 
quency break at Vb = 33.5^3 q mHz and a QPO with 
centroid frequency vqpo = 113.2 ±1.5 mHz in the power 
spectrum of M82 X-l. The QPO frequency is similar to 
that found from the RXTE observation of 1997 July 21, 
while it is about a factor of two larger than that found in 
the first XMM-Newton observation of 2001 May and an 
RXTE observation of 1997 February 24 (SM03). Thus 
the QPO frequency of M82 X-l is variable. 

The M82 X-l PDS is flat below the break frequency 
and then falls off above the break as a power law with 
a prominent QPO on top of it. This form of the PDS is 
characteristic of XRBs in their high/soft state or inter- 
mediate state (e.g., McClintock & Remillard 2003; see 
also Fig. H for a comparison). PDSs of many BH XRBs 
display low frequency QPOs in the frequency range of 
~ 0.05 - 30 Hz (Remillard et al. 2002). The QPO fre- 
quency is strongly correlated with the break frequency 
(Wijnands & van der Klis 1999). The break frequency is 
always lower than the QPO frequency by at least a factor 
of three (see, Wijnands & van der Klis 1999). M82 X-l 
follows this correlation as its QPO frequency is about a 
factor three larger than the break frequency. Thus the 
PDS of the ULX is very similar to the XRB PDSs in 
their high/soft or transition state. 

Since the characteristic time scales of accretion pow- 
ered sources scale with BH mass, the similarity of the 
PDSs of M82 X-l and XRBs provides the opportunity 
to determine the ratio of their BH masses. However, as 
stated above, the XRBs show a large range in their low 
frequency QPO and break frequencies. Therefore, the 
frequency scaling factor between any two power spec- 
tra of two XRBs will lead to an incorrect determina- 
tion of the ratio of their BH masses. Fortunately, the 
PDS features (break and QPO frequencies) are well cor- 
related with the energy spectral parameters. The low 
frequency QPOs are related to the flux of the PL en- 
ergy spectral component. The variable frequency QPOs 
in XTE J1550-564 and GRO J1655-40 appear only when 
the PL contributes more than 20% of the 2—10 keV 
flux (Sobczak et al. 2000a and references therein). 
GRS 1915+105 also shows similar behavior (Muno et al. 



1999). The PL index is well correlated with the centroid 
of the low frequency QPOs in GRS 1915+105 (Vignarca 
et al. 2003) and XTE J1550+564 (Sobczak et al. 2000a, 
2000b). The photon index increases with the QPO fre- 
quency until it saturates. The frequency scaling factor 
between the low frequency QPOs of two XRBs with sim- 
ilar PL photon indices can be used to determine the ratio 
of their BH masses (TF04; FT04). Thus it is crucial to 
determine the shape of X-ray spectrum and detect both 
the QPO and the low frequency break in order to make 
sure that the detected QPO lies at a low frequency, but 
above the break frequency, and on top of a PDS power- 
law continuum. 

The long XMM-Newton observation of M82 X-l has 
yielded all 3 pieces of information required to deter- 
mine its BH mass. The photon index of the 3 — 12 keV 
PL is 2.0 ± 0.1 as determined from the high signal-to- 
noise PN data. XRBs show a range of PL photon in- 
dices from r ~ 1.6 in their low/hard state to T ~ 2.5 
in their high/soft state. Thus M82 X-l is likely in 
a state intermediate between the high and low states. 
Sobczak et al. (2000a) studied the spectral behavior 
of XTE J1550-564 during its 1998-1999 outburst us- 
ing 209 pointed observations with RXTE. They used a 
model consisting of a multicolor disk blackbody and a 
PL to fit the spectra and found that the spectra grad- 
ually steepen from T ~ 1.5 to - 2.9 for XTE J1550- 
564. The low frequency QPO is also found to vary from 

0. 08 - 13 Hz. At a PL photon index of 1.91 (2.12), 
XTE J1550-564 shows a low frequency QPO at 0.81 Hz 
(1.6 Hz) (Sobczak et al. 2000a; 2000b). Noting that 
the BH mass scales inversely in proportion to the QPO 
frequency and the BH mass of XTE J1550-564 is about 
10M Q (Titarchuk & Shrader 2002, Shrader k, Titarchuk 
2003), we can determine the BH mass of the M82 X-l 
as M BH (ULX) ~ ^ QPO (XTE J1550-564)/i/ QPO (ULX) 
xM BH (XTE J1550-564) - 70- 150M©. Another XRB 
GRS 1915+105 shows a low frequency QPO with its cen- 
troid at 0.642 + 0.004 Hz (1.609 + 0.007 Hz) when its PL 
photon index is 1.88 ± 0.04 (2.16 ± 0.04) (Vignarca et al. 
2003). Noting that that the BH mass of GRS 1915+105 
is 14 ± 4M© estimated dynamically from the IR spec- 
troscopic observation of its companion star (Greiner et 
al. 2001), the BH mass for M82 X-l is M BH (ULX) 

- 50 - 260A/O. Thus, scaling the mass of M 82 X-l 
from the T — vqpo relation calibrated by the dynami- 
cal masses of several galactic BHs, we find Mbh (ULX) 

- 50 - 260M©. 

The QPO frequencies and photon indices for M82 X- 

1, XTE J1550-564 and GRS 1915+105 correspond to a 
region of steep dependence in the T — vqpo relation (see 
FT04). It is therefore likely that the uncertainty in the 
calibration of the T — vqpo relation may introduce an ad- 
ditional factor of two uncertainty in the determination of 
the BH mass of M 82 X-l. Thus the true BH mass of 
M 82 X-l is most likely in the range of 25 - 520M©. 
This is at least a factor of 2 lower than the previous es- 
timate by FT04, which can be explained by the lower 
values of the QPO frequencies for M82 X-l in our pa- 
per compared with those in FT04. The above estimate 
of the BH mass is based on the assumption that M 82 
X-l also follows the T — vqpo relation found for XRBs. 
XMM-Newton and RXTE observations of M 82 X-l have 
revealed that the QPO frequency is variable (SM03). The 
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photon index is also variable (FT04) but appears to be 
inconsistent with the T — vqpo relation for XRBs. How- 
ever, there may be systematic errors in the photon indices 
measured with XMM-Newton and RXTE due to contam- 
ination from nearby sources, as indicated by large appar- 
ent changes in the effective absorption column. Chandra 
ACIS observations of M 82 X-l were used to suggest 
an effective column in the range (5 — 9) x 10 20 cm~ 2 
(SM03), while XMM-Newton and RXTE observations 
have been used to suggest large columns in the range 
2 x 10 21 - 8 x 10 22 cm" 2 (see Section 3.2, FT04). Both 
XMM-Newton and RXTE observations suffer from con- 
tamination by nearby sources. A small contribution from 
nearby sources may significantly alter the spectral shape 
of the ULX and a varying contribution may lead to the 
observed changes in the absorption column. Thus, based 
on the currently available X-ray data on M 82 X-f, it 
is neither possible to prove nor rule out the Y — vqpo 
relation similar to that found for XRBs. 

M82 X-f is one of the brightest ULXs. During the 
second XMM-Newton observation, its 2 — 1 keV lumi- 
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Fig. 1 — PDS of the ULX M82 X-l and an XRB XTE J1550-564 having very similar PL photon indices (r ~ 2.0). Left: PDS of M82 
X-l derived from the EPIC PN and MOS data above 2 keV. The white noise level expected (~ 1.4) from the poisson errors has not been 
subtracted. The frequency resolution is 7.8 mHz. The best-fitting model comprising a broken power-law and a Lorentzian is shown as 
a thick line. Right: PDS of the XRB XTE J1550-564 derived from the RXTE observations of 1998 September 10. The best-fit model 
consisting of a broken power law and two Lorentzians is shown as a thick line. 
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Fig. 2. — (a) Ratio of EPIC PN data and the best-fit absorbed PL model showing the iron Kq line, (b) The PN spectrum and the 
best-fit absorbed PL+Gaussian line model (upper panel) and their ratio (lower panel). 



